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ABSTRACT 

The   elastic    differential    scattering    cross    sections    for 

+         +  +         +  o 

H    ,    H      +   Ar    -*   H    ,    H      +  Ar    were   measured   for    angles    of    39      to 

50      with    incident    lab    energy    of    100   eV .      The   apparatus    utilized 

the  properties    of   axially   symmetric,    non-uniform  magnetic    fields 

(large    thin    solenoid)    to   focus    all    scattered    ions    with    vector 

momentum  p    within    Ap    onto   a    detector    placed    down    the   magnetic 

axis.      Thus,    the   solid   angle   was    increased   by   a    factor    of   about 

300    over    conventional    methods    with    comparable   resolution    in    A0 

1   Q  Q 

(about    1    ).       A   typical    value   for    CT       (43    )     is    8.2   x    10  cm    . 

H+ 
The   H      data    was    compared  with    data    of   Magnuson,    et.al.,    via   a 

scaling    law   and  was    in    excellent    agreement.      A    typical    value 

,     o  ,  -20         2 

for    (J       (43    )     is    2.4   x    10  cm    .       A   simple   refinement    of   a 

H2 

+ 
classical    model    for    H      scattering   by    Feist    relates    these   cross 

sections    by  <T       (Q)    =   P  0"       (0)    ~   — % —  Lcr     ,(©)]      where   P    is    an 

Hj  H+  TTi>2E  H 

^  o 

orientation    factor,    D   the   binding    energy,    and  I    the    internuclear 
separation.       This    expression    relates    the   cross    sections    to   within 
the   accuracy    of    the   experiment    thus    verifying    the   classical 
appr  oximations . 
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I.   INTRODUCTION 

Atomic  processes  such  as  the  scattering  of  complex  ions 
by  heavy  atoms  are  properly  described  by  quantum  mechanics. 
Unfortunately,  for  a  many -body  problem  the  mathematics  becomes 
untractable.   In  order  to  produce  useful  solutions  approximations 
must  be  introduced.   One  class  of  approximations  involves  the 
use  of  semi -classical  models  of  the  collision  process.   By 
sacrificing  rigor,  useful  solutions  are  obtained  that  still 
provide  insight  into  the  physical  process  involved. 

The  classical  model  investigated  here  is  used  to  calculate 

+  . 
the  differential  scattering  cross  section  of  the  H   ion  by  a 

spherically  symmetric  potential.   The  theory  was  developed  by 
Feist  [l]  based  on  the  earlier  work  of  Bates,  Cook,  and 
Smith  [2]  as  follows.   First,  calculate  the  differential 
scattering  cross  section  of  a  single  proton  from  a  spherically 
symmetric  potential.   Then  calculate  the  probability  that  the 
second  proton  is  oriented  with  respect  to  the  first  such  that 
the  H   molecule  is  not  disrupted  by  the  collision.   The  dif- 
ferential scattering  cross  section  for  H   is  the  product  of 
this  orientation  probability  times  the  cross  section  for  the 
single  proton. 

The  purpose  of  the  experiment  reported  here  was  to  test 
the  theory  as  developed  by  Feist.   This  was  done  by  measuring 
the  differential  scattering  cross  section  of  H   by  Argon  and 


W      by    Argon   at    incident    energies    of    100    eV   and    laboratory   angles 
from    39      to   49    .      Thus,    the    orientation   probability    was    deter- 
mined  experimentally    and   compared  with    Feist's    calculation. 
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II.       THEORY 

A.       THE   MODEL 

This    classical    model    is    applicable    to    the   scattering   of    the 

+ 
H      ion    over    a    range    of    incident    energies    E      (in    laboratory 

coordinates)    by    a    spherically    symmetric   potential    such   as    heavy 

noble   gas    atoms.       The   reaction    studied  was    H„    +   Ar    **   H _    +   Ar 

where    the   underline    indicates    the    fast   particle    in    laboratory 

coordinates.      The    two  protons    of    the  H_    ion   have   an    equilibrium 

separation    in    the    ground    state    of    1.06    angstroms    and   a    disasso- 

ciation    energy    of   2.648    eV .      The    lower    energy    limit   where    this 

model    could   no    longer    be   expected    to   apply   would   be   where    the 

incident    ion    energy    becomes    comparable   with    the    ion    binding 

energy.      This    requires    that    E    >>  2.648    eV   or    E     "7    30    eV.      The 
^  M  o  o 

upper    energy    limit   where    this    model    could  no   longer    be   expected 
to   apply    is    governed   by    the   uncertainty   principle.       In    this 

application    it    is    required    that    the    de   Broglie  wavelength   be 

>  h 

much    less    than    the    distance    of   closest   approach    or    r   9?  K ^   =   

*^  o    D    mv 

o 

which  may  be  written  as  r   >  .   Now  r   may  be  further 

I  2mEQ 

related    to   E      by    considering   a    head-on    collision   with    the   shielded 
o 

coulomb  potential    for    Argon    (this    potential    is    discussed    in    Section 

II    B,    and   Appendix   A).       Considering    only    the    dominant    term   of   the 

potential    one    can   write 

2         -G!r 

E      -    Ae 

o  r 

o 

where  A   and  0L   are  parameters    of    the  potential.       Substituting 
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values  for  the  constants  one  obtains 

Eo  =  14.6  -\- 
o  2 

r 
o 


-1.68  r 

and   E   =  4ll  —  e 
o       r 
o 


with    E      in   units    of    eV   and   r       in   units    of   Bohr    radii.      The    value 
o  o 

A 
of   E      that    satisfies    these    equations    is    E      =    1.13   x    10      eV .       This 
o  o 

shows    that   a    classical   model    for    the    scattering    of   H      should   be 

valid    for    incident    energies    up    to   about    10  KeV . 

+ 
In   a    scattering    event    the  H      ion   may    be    disrupted   due   to    the 

following    two   effects.       First,    consider    the   two  protons    approach- 
ing a    spherically    symmetric   potential   with    the   same    impact   para- 
meter,   b.       Using    the   concept    of    an    impulse   approximation,    con- 
sider   each   proton-target   collision    separately,    yielding   a    de- 
scription   in    terms    of    two    two-body    collisions.      Due   to    the 
spherical    symmetry    each  proton    scattered    is    confined    to   a   plane 
defined   by    the   scattering  axis    and   a    line    from    the   scattering 
center    to    the  proton.      The   angular    difference   between    the 
scattering   planes,    A0   will    result    in   a    difference    in    velocities 

of    the    two   protons    after    scattering    such    that    their    relative 

+    . 
energy   may    be   greater    than    the   binding    energy    and   the   H?    ion 

will    be    disrupted.       Figure   1    illustrates    this   process. 

Second,    consider    the  protons    oriented   such    that    they   will    be 

scattered   in    the    same  plane   but   with    different    impact   parameters. 
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This    will    result    in    a    difference    in    scattering  angle   AO   and 
therefore   again   a    relative   velocity    between    the    two  protons. 
The  relative   energy    may   again    exceed   the    binding    energy    so    that 
the    ion   will    be   disrupted.       Figure   2    illustrates    this    process. 

In    the   actual    case    these   two  processes    are   not    separable 
and   disruption    may    occur    due   to   the   combined   effect    of   A0   and 

A9.      The   necessary    condition    for    the    ion    remaining    intact    after 

2 
collision    is    D>%/U(AV)       where   D    is    the   disassociat ion    energy, 

u,    =  M   /2    the   reduced  mass    of    the   two  protons,    and   AV    is    the 
p 

difference    in    velocities    due   to    both    A0   and    AO.      Whether    the 
H9    +   Ar    -*   H      +   Ar    reaction   will    occur    is    then    a    function    of   the 
orientation    of    the    two  protons    before    (or    equivalently    after) 
the   collision.      The   differential    scattering    cross    section    for 


the  H      ion,    <J      (0)    may    be   calculated  as    follows.       For    E  >>  D 

H2 
treat    the   collision    of    the  H_    ion   as    two    events    consisting    of 

the   scattering    of  proton    "A"   and  proton    "B"    (impulse   approxi- 
mation).      Calculate    the    differential    cross    section    for    proton 

"A",    CT       (9).       Calculate    the  probability,    P,    that   proton    "B"    is 

H+  2 

oriented  with   respect    to  proton    "A"    such    that    D>5g/i(AV)     .       It 


then    follows    that 


a  +(0)   =  p  a  +(0),  (l) 


h-  h 

The  electron  is  not  treated  as  a  scattered  particle  in  this 
model  but  may  be  thought  of  as  simply  providing  the  binding 
energy . 
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Bo   CALCULATION  OF  (J   (0) 

H 

The  calculation  of  O      (0)  was  done  on  an  IBM  36O-67  digital 

H+ 

computer    using   a   program   written    by    Feist    L  2 J  .      The  program 
utilizes    standard   numerical    techniques    to   solve    the   classical 
differential    scattering    equations.      The  potentials    used  are 
shielded   coulomb  potentials    from    Byatt    [3]    of    the    form 

2  -b   r 

V(r)    =  ^f-   (S    C   e      "    ) 

l         n      IJ 

with  EC      =1    but   C      and   b      otherwise  adjustable.       Byatt    lists 
n      n  n  n 

two  potentials    for    Argon,    one  with   and   one   without    correction 

for    polarization.      The    two  potentials    yielded   cross    sections 

+ 
for    H      that    were  within    13%   of   each    other    for    an    incident    energy 

of    100    eV   and   scattering   angles    from    38      to    50       (laboratory    co- 
ordinates).      Due    to    the    close   agreement    of    the    two  potentials, 
only    the   one    corrected   for    polarization    is    presented  here. 


By  a 


tt  indicates  the  potentials  should  be  accurate  to  within  5% 


■b  r 


for   £    C      e  >.01.       For    the    region   where   the  potentials    were 

n    n  "V 

used  it  was  noted  that  S  C_  e    >.052.   Thus  it  was  felt  that 


n 


n 


these  potentials    should   be   satisfactory    for    calculating  (J       (©) 

H 


C.       CALCULATION  OF   THE   ORIENTATION   PROBABILITY 

+  ... 

The   energy    criterion    for    the  H_    ion    remaining    intact   may    be 


put    as 


D  =  £    (AV 

2    v       max 


,         M 

2   =  J2L 


7-    (AV  ) 

k      v       max' 


or 


AV 


max 


4d 

M 


(2) 
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Considering  the  case  where  the  target  mass  is  much  greater  than 

+  + 

the  mass  of  the  H   ion,  the  H„  ion  has  energy  after  collision 

2 
E   =  M  V    both  in  the  laboratory  and  center  of  mass  systems, 
o     p  o 

Substituting  this  into  equation  (2)  gives 


AV 


max 


4d 

E 


V 


Now  it  was  shown  in  Section  II  A,  that  AV     =  f(b,  Ab,  A0). 

max  v  ' 

Therefore,    for    a    fixed    incident    energy   and    impact  parameter, 
there   exists    a    function    g    such    that    A0   =    g(Ab).       This    defines 
a    closed   curve    in   a   plane  perpendicular    to    the   scattering   axis 
that   proton    "B"    must   pass    through    in    order    for    the    ion    to   remain 
intact.      The   area    enclosed   by    the    A0    =   g(Ab)    curve    is    designated 

by    S. 

+ 
The   axis    of    the  H„    ion    is    randomly    oriented   with  proton    "B" 

having    no  preferred  position   with   respect    to  proton    "A".      The 
location    of  proton    "B"   will    therefore   define   a    sphere   about  pro- 
ton   "A"    of   radius   i  ,    the    internuclear    separation.      The    orientation 
probability    is    then    given    by    the   ratio   of   S    to    one-half    the   area 
of    this    sphere 


P    = 


2rrl 


2     ' 


(3) 


The  justification  for  using  the  area  of  a  hemisphere  in  the 
denominator  of  equation  (3)  is  that  the  "other  half"  is  re- 
dundant.  No  exact  solution  has  been  found  for  equation  (3)  and 
the  following  approximations  were  made. 
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First,  it  has  been  assumed  that  the  target  is  infinitely 


4 


heavy    so    that    the   speed   of    the  H      ion    is    unchanged   by  the    col- 
lision.      For    the   reaction    studied   experimentally,    this    is    valid 

in    that    M        =    2.032    amu   and  M         =    39-948    amu .       For    example,    a 

H*  Ar 

£  o 

100    eV   H      ion    scattered  at    45      has    an    initial    speed    of    9 • 80   x    10 

4 
m/sec    and   a    final    speed   of   9. 75   x    10      m/sec,       It    is    further 

assumed    that    the    two  protons    will    both    have    this    same   speed   after 
collision    and    that    disruption   will    occur    due   to   difference    in 
direction    only.       Feist   LlJ    showed   that    the   disruptive    effect 
of    the    difference    in    the   magnitudes    of   velocities    of    the    two 
protons    is    much    less    than    the    disruptive    effect    of    their    dif- 
ference   in    direction.       For    the    case  where    A0    =   0,    and   assuming 

A0    is    small,    AV  =   V   L  A0  I    rt         _..       Using    equation    (2)    this 

•"-"*'  o        max    A0   =    0  ^  v 


max 


becomes 


^■AGmax]A0    =    0 


4D 

E 
o  . 


(4 


Figure    3a    illustrates    this    case.       For    the   case  where    A©   =   0 

(or    equivalently    Ab   =   0)    and   assuming    A0    is    small, 

AV  =   V      sin    6    [A0         ] /.  _.       Using    equation    (2)    this 

max  o  max    Ab    =   0 


becomes 


£A0         Lk     •   n 
max    Ab    -   0 


4d 

E 


sin    0 


(5) 


Figure    3b   illustrates    this    case.       Using    the   classical    definition 

for    differential    scattering    I  2nbdbl     =    |  0"      (9)    2n    sin    0    d0|,     it 

H 


is    seen    that 


Ab 


Or       (0)    sin    0 
H+ 


m 


ax1 A0   =   0 


[AQmax]  A0    =   0    ' 


vTl- 


B   I  "  v0 


AV 


Proton  "B" 


max 


[A*max] 


A<£=0 


b 


Proton  HA" 


— 0 — 

Target  Atom 


Scatttring 
Axis 


Figure  3a.  Proton  velocity  difference  due  to  difference  in  impact 
parameter. 


v; 


*VQ   sin0 


max 


[**max] 


Ab'O 
Targtt  Atom 

Scattering  axis  into  paper 

Figure  3.  b.  Proton  velocity  difference  due  to  difference  in 
scattering  plane. 
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which,  upon  combining  with  equation  (4)  becomes 

a   (0)  sin  0 


Ab 


H 


max  A0  =  0 


4D 
E 

L  O 


(6) 


Equations    (.5)    and    (6)    give   maximum   values    for    either    Ab   or    A0 

+ 
such    that    the  H_    ion    is    not    disrupted.      The    A0   =   g(Ab)    curve 

must    therefore  be    enclosed   by    these   maximum   values.      An    upper 
bound   on    the  area,    S,    enclosed   by    the   A0    =   g(Ab)    curve    is    there- 
for e 

S  =    [2(Ab  )Anl    _    J[2b(A0  )  ^         J. 

max  x      max'A0    =   0  v       max'Ab   =   0 

Substituting    equations    (5)    and    (6)    into    the   above   and  using 
equation    (3)    an    upper    bound    on    the   orientation   probability    is 
found   to   be 


max 


«  CT       (e) 
J2-         H+ 

E  o 

o  .  2 


(7) 


and 


*    +(0) 


H- 


max        ^2E  H+ 


Due    to    the    simple   form    of   equation    (8),    it    is    calculated   numeri- 
cally   in    the   same   computer    program   used    to    calculate  0"       (©).       It 

H+ 
should   be    further    noted   that    D    is    the    disassociation    energy    for 

ground   state   and  I    the   corresponding    equilibrium   distance.      Ex- 

+ 
cited   electronic    states    for    the   H      molecule   are   either    unbound 

or    have   very    shallow   minima    in    their    potential    curves.       Therefore, 

if    the   H      ion    is    excited   to   higher    electronic    states    during    the 

collision   process    it    is    almost    certain    to   be   disrupted    due    to 
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velocity  differences  after  the  collision.   If  the  H9  ion  is 

excited  to  higher  vibrational  levels,  then  D  would  be  decreased 

and  I    increased.   The  molecule  could  still  remain  intact  but 

would  have  a  smaller  value  of  P  and  hence  CT   (©). 

H2 

D.       COMPUTER   CALCULATION  OF  CJ       (0)    AND  a       (0) 

H+  H+2 

Figure    4    displays    the   results    of    the    computer    calculations 

for    the   differential    scattering    cross    sections.      One   set   of 

+  +  ... 

curves    is    for    the   H_   +   Ar    ~*   H      +   Ar    reaction   at    incident    energies 

of    100,300,    and   500    electron-volts.      The   second   set    of    curves 

+  +  . 

is    for    the   H      +   Ar    -»   H      +   Ar    reaction   also   at    incident    energies 

of    100,    300,    and   500    electron-volts.      The   H      curves    were    com-    . 
puted  using    equation    (8)    and  represent   a   predicted   upper    limit 
on    the   scattering   cross    section. 

It  is  noted  that  all  of  the  curves  are  monotonically  de- 
creasing functions  of  scattering  angle  and  reflect  no  atomic 
str uctur  e . 
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Figure  4.  Log  of  scattering  cross  section  vs.  scattering  angle  for 


H*+  Ar 


H*  +  Arand  H2*  +  Ar 


H2*  +  Ar. 
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III.   EXPERIMENTAL  APPARATUS 

A.   GENERAL  DESCRIPTION. 

The  experimental  apparatus  used  to  measure  the  cross  sections 

+    +  +    + 

for  the  H  ,  H   +  Ar  -  H  ,  H   +  Ar  reactions  is  discussed  in  detail 

by  Bush  L 4] .   An  overall  schematic  of  the  apparatus  is  shown  in 

Fig.  5-      The  operation  of  the  apparatus  was  as  follows.   Three 

+    +    + 
species  of  hydrogen  ions  (H  ,  H  ,  H  )  were  produced  in  the  duo- 

plasmatron  ion  source  and  were  accelerated  electrostatically 
to  the  desired  energy.   The  ions  entered  a  mass  analyzer  where 
the  desired  ion  species  was  bent  through  a  45   angle  to  travel 
down  the  geometric  axis  of  the  apparatus.   The  selected  beam  was 
focused  and  entered  the  scattering  cell  where  the  target  gas, 
in  this  case  Argon,  was  contained.   The  scattering  cell  was 
centered  in  a  large  focusing  magnet  that  produced  a  known  non- 
uniform but  axially  symmetric  magnetic  field.   The  trajectory  of 
an  ion  with  a  given  vector  momentum  could  then  be  calculated  and 
the  detector  was  appropriately  placed  to  collect  this  ion.   Due 
to  the  finite  detector  size,  ions  scattered  within  a  finite  A© 
(on  the  order  of  1  )  were  also  detected.   It  is  further  noted 
that  essentially  all  the  ions  scattered  at  an  angle  0  were  col- 
lected regardless  of  their  azimuthal  angle  of  departure.   Cross 
section  measurements  at  different  scattering  angles  could  be  made 
in  two  ways.   Either  the  detector  could  be  fixed  on  the  magnetic/ 
geometric  axis  and  the  focusing  magnet  current  varied  to  select 
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the  desired  scattering  angles,  or  the  focusing  magnet  current 
could  be  fixed  and  the  detector  moved  along  the  magnetic/geometric 
axis  to  select  the  desired  scattering  angles.   In  this  experiment, 
the  latter  method  was  used.   The  detector  system  consisted  of  a 
magnetic  electron  multiplier  whose  output  was  fed  into  an  elec- 
trometer.  The  entire  system  was  maintained  at  a  vacuum  on  the 
order  of  10    torr  during  operation  by  two  diffusion  pumps. 

B.   THE  DUOPLASMATRON  ION  SOURCE. 

The    duoplasmatr on    ion    source   was    constructed   and    described 
in    detail    by    Carter    L  5J  ■•       Figures    6    and   7   display    the   construction 
of    the   source   and   the   circuit    diagram    for    its    operation.       The 
operation    of    the    source   basically    consisted   of   maintaining    the 
plasma    discharge    in    the   Z-electrode    cavity   and    extracting    ions 
at    the    desired   energy    by   placing    this   potential    on    the   anode. 
Focusing    the    extracted   beam  was    accomplished  by    an    electro-magnet 
surrounding    the   Z-electrode   cavity    and  a    conventional    Einzel    lens 
and   deflector    stack. 

To    increase    the   beam    intensity    an    acceleration-deceleration 
system  was    utilized.      This    consisted   of   extending    the    extraction 
lens,    L    ,    to    enclose    the   entire    focusing    stack.      This    enclosure 
was    called    the   "sock".       (See   Fig.    7.)       By    floating    the    sock  at    a 
negative  potential    V    ( typically-500    volts)    and   by    increasing    the 
potentials    of    the    interior    lenses    and   deflectors    by    the   same 
amount,    the   extracted    ions    transversed   this    region   with    an    energy 
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(E  +  V)e,  yielding  a  higher  space  charge  limit  for  the  beam 
density.   Upon  exiting  from  the  sock,  the  hydrogen  ions  were 
decelerated  to  energy  Ee.   This  arrangement  increased  the  beam 
intensity  at  the  scattering  cell  by  a  factor  of  about  40 . 

A  typical  maximum  beam  intensity  for  mass  analyzed  protons 

of  100  eV  delivered  to  the  scattering  cell  with  0,5  cm  entrance 

-7 
aperture  at  80  cm  from  the  source  was  1  x  10    amperes.   For 

+ 
H   ions  under  the  same  conditions  the  maximum  beam  intensity  was 

-7 
on  the  order  of  2  x  10    amperes,   The  energy  spread  of  the  beam 

was  estimated  to  be  about  2  eV  for  a  100  eV  beam. 


C.   THE  MASS  SPECTROMETER. 

The   mass    spectrometer    used    to   select    the   desired    ion   was 
constructed,    calibrated  and   described    in    detail    by    Strohsahl    \_6] 
Figure   8    is    a    schematic    of    its    arrangement    in    the   overall    appa- 
ratus.      The   mass    spectrometer    utilized  an    electromagnet    to  pro- 
duce  a    magnetic    field    of   up    to    4    kilogauss    confined   between    two 
40      wedge-shaped  pole  pieces.       This    was    sufficient    to   mass 
analyze   atomic   masses    up    to   40   with    energies    up    to   1   KeV.      The 


resolution    (- — )    was    about    40 .      A    bending   angle   of   45      was 

v  Am 

utilized  and  this  altered  the  selected  ion  beam  trajectory  from 
the  axis  of  the  duoplasmatr on  to  the  geometric  axis  of  the  main 
scattering  apparatus.   Operation  of  the  mass  analyzer  consisted 
of  merely  of  setting  the  correct  magnet  current  to  focus  the 
desired  ion  into  the  scattering  cell.   The  pole  faces  and  region 
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where  undesired  species  were  "dumped"  were  coated  with  aquadag 
to  prevent  static  charges  from  building  up.   Even  so,  this  region 
required  thorough  cleaning  after  about  20  hours  of  duoplasmatr on- 
analyzer  operations. 

D.   THE  SCATTERING  CELL. 

The  scattering  cell  was  designed  by  Bush  L 4] •   Figure  9  shows 
a  schematic  of  its  construction.   The  functions  of  the  scattering 
cell  are  to  (1)  allow  the  collimated  incident  ion  beam  to  pass 
through  the  scattering  region  and  be  measured  by  the  beam  collector, 
(2)  contain  the  target  gas  in  the  scattering  region  and  measure 
its  pressure,  and  (3)  allow  scattered  ions  to  escape  over  as  large 
a  range  of  scattering  angle  and  azimuthal  angle  as  possible.   The 
design  of  the  scattering  cell  was  a  compromise  between  these  re- 
quir  ements . 

The  incident  beam  enters  the  scattering  cell  through  an 
adjustable  aperture  (diameter  from  0.1  to  1.0  cm)  which  was 
usually  set  at  0.5  cm.   The  beam  passed  through  the  front  section 
which  contains  the  scattering  region.   Now  the  great  majority  of 
the  incident  ions  will  not  be  scattered  at  all  and  will  continue 
on  into  the  rear  section  where  they  are  collected  either  by  the 
rear  section  wall  or  the  beam  collector.   The  incident  beam  was 
measured  by  two  Keithly  Model  4l0  Micr o-Micr oammeter s .   The  beam 
collector,  which  collects  about  90%  of  the  beam,  had  three  grids 
that  prevented  slow  ions  from  being  collected  and  suppressed 
secondary  electron  emission.   The  beam  collector  assembly  could 

be  swung  off  axis  allowing  a  clear  path  through  the  scattering  cell 

for  system  alignment  and  detector  calibration. 
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The  target  gas  was  bled  into  the  scattering  region  through 

a  Granville-Phillips  variable  leak  valve.   The  pressure  was 

measured  by  a  VG-1A  ion  gauge  manufactured  by  Vacuum  Industries. 

Since  operation  of  the  ion  gauge  interfered  with  the  incident 

beam  measurement;  the  pressure  was  measured  before  and  after 

each  run  and  averaged.   Vacuum  Industries  lists  the  correction 

factor  for  Argon  (relative  to  dry  air)  as  0.617.   That  is, 

P      ■  =  0.617  x  P.  ...      . 
actual  indicated 

The  rear  section  of  the  scattering  cell  was  adjustable  along 
the  axis  of  the  apparatus  and  thus  regulated  the  target  thickness 
between  0.1  and  0.8  cm.   This  was  set  at  0,3  cm  for  the  present 
expei iment . 

The  scattering  angle  limits  were  defined  by  the  slope  of  the 
front  and  rear  section  shoulders  of  49   and  3°   respectively. 
This  range  of  scattering  angles  could  be  extended  somewhat  by 
considering  an  effective  target  thickness.   Bush  L4J  describes 
this  process.   Ideally,  ions  scattered  into  all  36O   of  azimuthal 
angle  could  be  collected  at  the  detector.   However,  interference 
by  mechanical  support  required  for  the  scattering  cell  reduced 
the  available  azimuthal  angle  to  290  . 

E.   THE  FOCUSING  MAGNET. 

The   focusing   magnet    consisted    of   five    concentric    spools    of 
copper    wire    in   a   parallel-series    arrangement.      Each    spool   was 
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water-cooled  by  a  \    inch  copper  tubing  jacket.   A  Sorenson 
Nobatron  DCR  150-15  power  supply  maintained  the  desired  current 
setting  for  periods  of  several  hours  with  essentially  no  drift. 

The  magnetic  field  on  the  center  axis  was  carefully  surveyed 
by  Kelly  L7J  using  a  Hall  probe.   Figure  10  shows  the  result  of 
these  measurements.   Kelly  also  showed  that  the  magnetic  field 
was  symmetric  about  the  coil  axis  to  within  0.1%.   This  allowed 
the  axis  of  the  apparatus  to  be  aligned  with  the  coil  axis  so 
that  the  magnetic  axis  coincided  with  the  geometric  apparatus 
axis  . 

In  order  to  better  understand  the  function  of  the  focusing 

magnet,  it  is  instructive  to  consider  the  simpler  case  of  a 

uniform  axially  symmetric  magnetic  field.   This  would  correspond 

to  the  entire  apparatus  being  enclosed  in  a  long  solenoid.   For 

an  ion  of  charge  q  with  momentum  p  and  scattered  at  an  angle  0 

the  trajectory  is  a  helix  such  that  after  one  orbital  revolution 

the  ion  will  again  intercept  the  axis  at  some  distance  Z  . 

^  o 

Using  these  boundary  conditions  the  Lorenz  force  equation 

F  =  q  v  x  B  can  be  integrated  to  give  Z   =  —  p  cos  0.   A  de- 

o    qB  r 

tector  placed  at  Z  will  then  detect  the  scattered  ion.   (See 

r  o  v 

Fig.  11.)   This  equation  also  shows  the  two  ways  of  detecting 
scattered  ions  over  a  range  of  0  and  p.   Either  maintain  B  con- 
stant and  vary  the  detector  position  or  keep  the  detector  fixed 
and  vary  B. 
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Now  in  order  to  measure  a  differential  scattering  cross 
section,  the  solid  angle  that  the  finite  size  detector  "sees" 
must  be  known.   For  the  axially  symmetric  case  this  is 
dCl    =    2tt  sin  0  d0.   To  evaluate  d9  consider  the  various  para- 
meters affecting  Z  ,  the  axis  intercept  distance.   A  functional 
relationship  exists  such  that: 

Zq   =    f(E,O,Bs/0,775C) 


wher  e 


E   =    energy    of    incident    ion 

0   =    angle    of    scatter 

B  =    magnetic    field 

p   -    radial    target    displacement    from    the    system 
axis 

77    =    displacement    of    target   atom   from    the 

scattering   center    along    the    system   axis 

C   =    constant    for    each   reaction    including    mass 
ratios,    inelastic    energy    losses,    etc. 


Differentiating    this    expression    gives 


A    Z      = 
o 


'dZoN 

~6E. 


AE    + 


s_5 


AG    + 


d_5 

dB 


'dz 


AB    + 


bp 


A/0     + 


For    a    specific    measurement    AB  and    AC   are   zero.       Since   analytic 

expressions    for    the   various    derivatives    are   not   known,     they   must 

be   evaluated    numerically    by    varying   E,    0,   p,    and  77,  .       Bush,    Cook, 

Heinz,    and   Rodeback   L  8  J    estimate    that    A/0    and    A??    which   are    due 

to    the    finite    volume    of    the   scattering    region    affect    AQ    by    no 

more    than    10%.       Approximating    AjD    and   AT?    as    zero    the   above   ex - 

/dz_\ 


pression  can  be  rearranged  to  A9  = 


'bQ 

az. 


AZ 


dE 


AE 
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The  effect  of  a  small  (2%)  energy  spread  was  studied  and  it  was 
concluded  that  the  "effective"  AO  was  not  changed.   The  angular 
resolution  of  the  apparatus  was  degraded  somewhat  in  that  not 
all  particles  scattered  within  AQ  of  0  were  collected  but  this 
was  compensated  by  the  fact  that  some  particles  scattered  outside 
of  A©  were  collected.   The  solid  angle  can  therefore  be  approxi- 
mated as  /.   \ 

dft    ~    2u    sin    0    l(l^)   |     AZQ 

AZ 
where   +_  — —  is    the   maximum    deviation    of    the    trajectory    inter- 
cept   from   Z  that   will    still    allow    the   particle    to    be    collected 
o  r 

by    the    detector.       AZ      can    be    calculated   by    considering    the   de- 
tector   geometry    and    the   angle   0      with   which    the  particle    inter- 
cepts   the   Z-axis.      As    an   example   of   how   the  approximate   expression 
for    the    solid   angle    is    applied    if   an   analytic    expression    for 

Z      =    f(0)     is    known,    consider    again    the   uniform   axially    symmetric 

2rr 
magnetic    field.       Recalling    that    Z      =   —   cos    0,    it    follows    that 

o         qB 

-    p  o 

For    the    non-uniform   axially    symmetric    field,    the  procedure 
for    computing    trajectories    and    solid   angles    was    essentially    the 
same   as    for    the   uniform    case   except    that    the    calculations    were 
done   numerically    on    an    IBM    360-67    digital    computer.       The   program 
that   accomplished    this    is    explained    in   Appendix  A.      A    very    brief 
outline   of    the  program    is    as    follows. 
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(1)  Express  the  magnetic  field  along  the  axis  as  a  12th 

13        n_1 

order    polinomial.       This    gives    B      =    (I)   £    C   Z  where 

z         v     '    i      n 

I    is    the   magnet    current. 

(2)  Using    this    polynomial,    a    vector    potential    is    found. 

(3)  Combining    the    vector   potential    with    the  Lorentz    force 
equation    gives    a    2nd    order    non-linear    differential 
equation.       This    is    expressed   as    two    simultaneous,    first- 
order    differential    equations    which   are    integrated   by 
the    DHPGC   Subroutine    of    the   Scientific    Subroutine 
Package   for    the    IBM    36O-67    computer.      The    trajectory, 

Z    ,    and   0      are   now   known. 
o  o 

(4)  The   above   procedure    is    done   over    a    selected    range   of 
scattering    angles    and   a   polynomial    expressing    0   =    f(Z    ) 

is    generated. 

(5)  From   0      and    the   detector    geometry,    AZ      is    calculated. 

v  '        o  o 

(6)  The  solid  angle  can  now  be  calculated  as 

is  present  due  to  the  construction  of  the  scattering 
cell. 
Figure  12  shows  a  series  of  trajectories  computed  for  H 
ions  with  incident  energy  of  100  eV  elastically  scattered  by 
Argon . 
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F.   THE  DETECTOR. 

The  detector  system  employed  in  this  experiment  consisted  of 
a  Bendix  Model  306  Magnetic  Electron  Multiplier  (MEM)  which  was 
fed  into  a  Keithly  Model  640  Electrometer,,   Figures  13,  14  and 
15  show  the  detector  geometry,  a  schematic  of  operation,  and 
a  circuit  diagram.   The  principle  of  operation  of  the  multiplier 
was  as  follows.   An  incident  ion  was  accelerated  to  about  1500  eV 
and  ejected  an  electron  from  the  tungsten  cathode  upon  impact. 
A  crossed  electric  and  magnetic  field  (alnico  permanent  magnets) 
imparted  a  cyclical  trajectory  to  the  electron(s)  down  the  dynode 
strip.   The  glass  dynode  and  field  strips  had  a  high  resistance 
coating  of  tin  oxide  and  antimony  that  sustained  the  electric 
field.   Each  time  the  electrons  struck  the  dynode  strip  multi- 
plications occurred  so  that  when  the  electrons  were  collected  by 

4 

the  anode,  a  gain  of  about  10   had  been  realized. 

Unfortunately  the  MEM  had  to  be  operated  in  the  focusing 
magnetic  field.   This  resulted  in  a  sideways  drift  of  the  elec- 
trons moving  down  the  dynode  strip  so  that  all  of  the  electrons 
generated  were  not  collected  at  the  anode.   This  in  fact  made 
the  MEM  gain  a  very  sensitive  function  of  its  position  in  the 
focusing  magnetic  field.   The  procedure  used  was  to  calibrate 
the  MEM  under  the  same  environment  that  existed  during  the  actual 

cross  section  measurements.   Specifically,  this  consisted  of: 

+ 
(1)  measuring  the  gam  using  H„  ions  with  a  focusing  magnet  cur- 
rent of  6.6  amperes  with  the  detector  position,  Z  ,    varied  from 
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30.5  to  47»5  centimeters  from  the  scattering  cell,  (2)  measuring 

+ 
the  gain  using  H   ions  with  a  focusing  magnet  current  of  4  °  8 

amperes  with  the  detector  position,  Z  ,  varied  from  28 .5  to  42.5 

centimeters  from  the  scattering  cell.   Since  the  focusing  magnetic 

field  was  a  known  function  of  current  and  position,  the  two  gains 

are  plotted  against  magnetic  field  in  Fig.  l6 „   The  gain  for  H_ 

ions  was  measured  before  and  after  the  cross  section  measurements 

and  was  found  to  remain  constant.   The  gain  for  H   ions  was 

measured  only  once. 

Due  to  the  steep  slopes  of  the  gain  curves,  a  small  error 

in  setting  detector  position  or  focusing  magnet  current  will 

result  in  a  large  error  in  the  detector  current  reading.   This, 

coupled  with  the  tedious  task  of  calibrating  the  detector  for 

every  species,  energy,  and  magnetic  current  that  a  cross  section 

measurement  is  to  be  made,  made  this  detector  system  less  than 

ideal.   For  further  experimental  work  with  this  apparatus,  a 

different  detector  system  is  being  considered.   This  may  consist 

of  counting  the  output  pulses  from  the  MEM  rather  than  integrating 

them  with  an  electrometer. 
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Figure  16.  Multiplier  gain  vs.  magnetic  field 
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IV.   EXPERIMENTAL  RESULTS 

A.   LABORATORY  CROSS  SECTION  MEASUREMENT. 

The  equation  used  for  calculating  a  differential  cross  section 

in  the  laboratory  was 

_    V^ 2 

d(T(0)         I  .x(yO   x   P         x   Corr  .Factor  )x    t   x   dO(0)    Cm    * 

where  I      =    detector    current   measured   by   an    electrometer 

G   =    MEM    gain 

I.  =  Incident  ion  beam.   Sum  of  the  currents  to 

1 

Beam  Collector  and  Rear  of  Scattering  cell. 

l  ft  ^ 

D    =    3 • 54  x  10    particles/cm  -torr  based  on  ideal 

gas  law. 

P         =    Pressure    in    scattering   cell    as    measured    by 
sc 

VG-1A    ion    gauge 
Corr .    Factor    =    .617;    corrects    ion    gauge   reading   for   Argon 
t    =    target    thickness.       Normally    equal    to    .3    cm 
in    this    experiment 
dO(0)    =    Solid   angle   computed   by    digital    computer. 
A   sample   calculation    for    the   H      +  Ar    reaction    showing    typical 
experimental    values    is 

4 


a       (Ll    Q°)    =    (1-85    x    10         amp.)/1.27   x    1CT 
+  ,      ,  -.  „-8  .  ,„    _,  ^1 


H'  (1.41    x    10""    amp.)(3.54   x    10      cm      -torr       ) 


(1.38    x    10"3torr    x    .6l7)(0.3    cm)(0.109    sterad) 

-17         2 
=    1.07   x    10       '    cm 
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B.   MEASUREMENT  OF  THE  H   +  Ar  -  H   +  Ar  CROSS  SECTION. 

The  differential  scattering  cross  section  for  this  reaction 
was  measured  from  39-2   to  50  .   Five  runs  were  made  and  the 
data  points  are  displayed  in  Fig.  17  along  with  the  computed 
cross  section  using  the  Byatt  potential.   Figure  l8  shows  the 
average  of  the  5  runs  along  with  the  computed  cross  section. 
Important  parameters  for  the  measurements  included: 

(1)  Incident  energy  of  H   ions  equal  to  100  eV. 

(2)  Focusing  magnet  current  equal  to  4.8  amps. 

(3)  Solid  angle  typically  equal  to  0.1  steradians. 

(4)  Argon  pressure  in  scattering  cell  typically  equal  to 

-3 
1  x  10  J  torr . 

Figure  l8  shows  that  the  experimental  curve,  while  having 
about  the  same  overall  slope  as  the  computed  curves,  differs  in 
magnitude  by  a  factor  of  about  6.   This  difference,  while  not 
catastrophic,  was  greater  than  expected.   First,  there  may  have 
been  a  systematic  error  in  the  experimental  apparatus.   A  pos- 
sible systematic  error  could  have  been  an  error  in  measuring  the 
MEM  gain.   The  gain  was  measured  with  an  incident  beam  several 
orders  of  magnitude  greater  than  the  beam  of  scattered  particles. 
Thus,  if  the  MEM  gain  were  a  (slowly)  increasing  function  of  input 
current,  the  measured  gain  would  be  greater  than  the  gain  for  the 
smaller  beam  of  scattered  particles.   This  would  result  in  the 
measured  cross  sections  being  too  small. 
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Figure  17.  Experimental  measurements  of  ^1*  +  Ar  -►  _H*  +  Ar  at 
100  eV. 
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Second,  the  potential  may  have  been  inadequate  for  this 
application.   Varying  the  potential  parameters  revealed  that  the 
computed  cross  sections  were  relatively  insensitive  to  small 
changes  in  the  potential  parameters.   For  example,  varying  the 
dominate  terms  in  the  potential  by  30%  such  that  the  collision 
was  made  softer  decreased  the  computed  cross  sections  by  a  factor 
of  only  about  2.   It  was  felt  that  the  potential  parameters  were 
at  least  this  accurate.   It  is  noted,  however,  that  the  potential 
is  Hartree-like  and  therefore  will  reflect  no  atomic  structure 
or  inelastic  process.   At  the  large  angles  where  this  experiment 

was  conducted,  the  collisions  are  relatively  hard  so  that  the 

+ 
Argon  atom  may  well  be  excited.   The  H   ion  would  then  have 

significantly  less  energy  and  momentum  after  the  collision  com- 
pared to  the  elastic  case.   Since  the  focusing  magnetic  field 
acts  as  a- momentum  analyzer,  only  elastically  scattered  H   ions 
were  collected  by  the  detector.   Thus,  the  presence  of  an 
inelastic  process  could  decrease  the  elastic  cross  section  by  a 
factor  of  6  and  not  be  accounted  for  in  the  Hartree-like  potential. 
The  work  of  Aberth  and  Lorents  L 11 J  tends  to  support  the 

second  analysis.   Aberth  and  Lorents  measured  the  scattered  signal 

+  o      , 

profile  of  He   on  Ar  at  10   and  hOO  eV  with  an  energy  analyzer. 

The  result  shows  that  for  He   on  Ar ,  an  inelastic  event  is  more 

likely  than  elastic  scattering  by  a  factor  of  about  2%.   The  energy 

lost  by  He   in  inelastic  events  was  on  the  order  of  10  to  15  eV. 
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Applying  this  to  the  case  of  H   on  Ar ,  it  seems  reasonable  to 
expect  protons  to  be  a  more  efficient  projectile  for  exciting 
Ar  than  He   and  therefore  could  account  for  the  factor  of  6. 
Furthermore,  inelast ically  scatter  protons  with  energy  after 
collision  of  85  -  90  eV  would  not  be  collected  by  the  detector. 

The  choice  between  the  computed  and  experimental  cross 
sections  was  accomplished  by  comparing  the  data  in  Fig.  l8  with 
data  supplied  by  Magnuson,  Carlston,  and  Snyder  [ 9] •   Magnuson 
et  al  also  measured  the  cross  section  of  H   on  Ar  although  their 
measurements  were  made  for  scattering  angles  less  than  20   and 
energies  ranging  from  150  eV  to  800  eV .   Meaningful  comparison 
of  the  data  was  possible  by  using  a  reduced  plot  outlined  by 
Smith,  Maichi,  Aberth,  Loients,  and  Heinz  LlOj.   This  reduced 
plot  consists  of  plotting  T  =  E©  versus  p  =  Qsin  0  (7(0).   By 
the  use  of  this  scaling  law,  cross  section  curves  obtained  at 
various  energies  are  reduced  to  a  single  line.   The  reason  is 
that  the  product  E0  represents  essentially  the  same  encounter 
with  the  atomic  potential.   By  comparing  the  data  in  this  manner 
it  was  found  that  the  two  experimental  curves  agreed  to  within 
a  few  percent  and  the  unaveraged  data  infact  overlapped.   Thus, 
it  was  concluded  that  the  experimentally  determined  cross  sections 
were  probably  correct.   It  is  further  noted  that  the  data  of 
Magnuson,  et.al.  covered  a  much  larger  range  of  T    =  E©  than  this 
experiment  and  does  reflect  atomic  structure  such  as  discussed 
in  the  preceeding  two  paragraphs. 
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C.   MEASUREMENT  OF  THE  H2  +  Ar  -  H_"  +  Ar  CROSS  SECTION. 

The  differential  scattering  cross  section  for  this  reaction 
was  measured  from  38-8   to  49-7  •   Five  runs  were  made  and  the 
data  points  are  displayed  in  Fig.  19=   The  average  of  the  5  runs 
is  displayed  in  Fig.  20.   Important  parameters  for  this  measure- 
ment included: 

(1)  Incident  energy  of  H   ion  beam  equal  to  100  eV . 

(2)  Focusing  magnet  current  equal  to  6.6  amperes. 

(3)  Solid  Angle  typically  equal  to  0.1  steradians. 

(4)  Argon  pressure  in  the  scattering  cell  typically  equal  to 

-3 
1  x  10  J  torr . 

From  Fig.  20  it  is  seen  that  the  overall  slope  of  the  experi- 
mental curve  is  essentially  flat.   That  the  cross  section  does 
not  decrease  significantly  with  increasing  angle  is  possibly  due 
to  the  following  instrumental  effect.   Some  particles  scattered 
outside  of  the  computed  solid  angle  may  reach  the  detector  through 
multiple  collisions  with  the  walls  of  the  apparatus  and  still  have 
sufficient  energy  to  overcome  the  detector  grid  bias.   The  sigi- 
nal  due  to  these  extraneous  particles  could  not  be  subtracted 
out  in  that  they  essentially  behaved  like  the  "real"  particles 
scattered  into  d  Cl.     That  is,  the  signal  due  to  extraneous 
particles  was  proportional  to  the  incident  H      beam  intensity  and 
to  the  pressure  in  the  scattering  cell.   Normally,  this  extraneous 
signal  is  negligable  compared  to  the  signal  from  particles  scat- 
tered within  d  f2.  However,  in  this  case  at  the  larger  angles  the 
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detector  signal  was  quite  small,  typically  on  the  order  of 

-12 
0.3  x  10     amperes.   Therefore,  perhaps  up  to  50  percent  of 

the  measured  signal  could  have  been  due  to  these  extraneous 

particles.   This  would  also  account  for  the  wider  dispersion  of 

data  at  the  larger  angles. 

D.   EXPERIMENTAL  RESULTS  AND  THE  CLASSICAL  MODEL. 

Since  it  could  not  be  assertained  with  certaintly  as  to 
whether  the  computed  or  experimental  cross  section  for  H   on 
Argon  was  more  accurate,  the  analysis  of  the  classical  model 
was  done  from  both  points  of  view. 

First,  it  was  assumed  that  the  computed  cross  section  for  H 
was  correct.   The  difference  between  the  experimental  and  com- 
puted cross  section  was  then  attributed  to  a  systematic  error 

+ 
and  the  experimental  H   data  was  normalized  to  fit  the  computed 

curves.   This  was  done  by  moving  the  H   data  up  by  a  factor  of  o. 

+ 
Now  the  H   on  Argon  cross  section  measurement  had  the  following 

+ 
similarities  with  the  H   measurement: 

(1)  Approximately  the  same  focusing  magnetic  field  (see  Fig.  16 ) 

(2)  Approximately  the  same  excursion  in  detector  position. 

(3)  Approximately  the  same  calculated  solid  angle. 

(4)  The  gain  curves  for  both  measurements  have  about  the  same 
shape  (see  Fig.  l6 ) . 

It  was  therefore  concluded  that  if  a  systematic  error  were  pre- 

+ 
sent,  it  should  have  affected  the  H   measurement  by  about  the 
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+ 
same  magnitude.   Therefore,  the  H_  data  was  also  increased  by  a 

factor  of  6„   The  results  of  this  normalization  are  shown  in 

+ 
Fig.  21 .   It  is  seen  that  the  H   on  Argon  experimental  data  pre- 
sented in  this  manner  lies  below  the  predicted  curve  by  a  factor 
of  about  l4« 

Second,  it  was  assumed  that  the  experimentally  measured  H 
on  Argon  cross  section  was  more  accurate  than  the  computed. 

Since  the  classical  model  predicts  that  the  cross  section  for 

+  +  n 

H   is  proportional  to  the  square  of  the  H   cross  section  (Eq.  o), 

+ 
the  predicted  curve  for  H   on  Argon  in  Fig.  21  would  be  moved 

,2      , 
down  by  a  factor  of  (6)   or  3°.   Figure  22  illustrates  this.   In 

this  case,  it  is  seen  that  the  experimental  measurement  for  H_ 

on  Argon  lies  below  the  cross  section  predicted  by  the  classical 

model  by  a  factor  of  about  2  to  3= 
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V.   CONCLUSIONS 

The  experimental  results  given  above  can  be  analyzed  in  two 
distinct  ways:  (1)  Assume  that  the  theoretical  H   +  Ar  curve  is 
correct  and  that  a  systematic  error  was  present  in  the  experi- 
mental measurements.  (2)  Assume  that  the  experimental  H   +  Ar 
cross  section  is  correct.   We  present  below  the  conclusions  based 
on  both  assumptions  and  our  reason  for  believing  that  method  (2) 
is  more  likely  to  be  correct. 

If  a  systematic  error  were  present  and  affected  the  measure- 

+       + 
ments  of  both  the  H   and  H   cross  sections  by  the  same  magnitude, 

then  it  was  seen  that  the  classical  model  predicts  a  cross  section 

for  H   too  large  by  a  factor  of  about  lq. .   That  the  actual  cross 

+ 
section  for  H   is  so  much  smaller  than  predicted  is  probably  due 

to  the  two  effects  previously  discussed.   That  is:  (1)  the 

orientation  probability  was  computed  by  considering  the  maximum 

area  enclosed  by  the  A0  =  g( Ab)  curve  and  (2)  the  model  neglected 

+ 
the  possibility  that  the  H   molecule  might  be  excited  during  the 

collision  process.   That  the  predicted  cross  section  was  so  much 

greater  than  observed  would  perhaps  make  the  classical  model 

+ 
useful  only  for  predicting  upper  limits  of  Hp  cross  sections. 

+ 
However  in  view  of  the  agreement  of  our  H   on  Argon  data  with 

Magnuson  et.al.  [  9*1  ,  this  is  considered  the  less  likely  of  the 

two  conclusions. 

If  both  experimental  measurements  are  correct,  then  it  is  seen 

that  the  classical  model  predicts  an  H      cross  section  too  large 
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by  a  factor  of  about  2  to  3.   This  would  make  the  classical  model 

+ 
useful  for  predicting  H   cross  sections  on  an  absolute  basis. 

This  is  considered  to  be  the  more  likely  of  the  two  conclusions. 

Assuming  that  the  classical  model  is  essentially  correct  the 

following  simple  refinement  is  proposed:   Recalling  equation  (3) 

it  was  seen  that  the  orientation  probability  could  be  written  as 

S 


P  = 


2v\i2 


where  S  was  the  area  enclosed  by  the  A0  =  g(Ab)  curve.   Since  S 
could  not  be  determined  exactly,  an  upper  limit  was  obtained  by 
considering  the  maximum  values  for  A0  and  Ab  respectively  when 
the  other  variable  was  zero.   This  in  fact  gave  four  points  on 
the  A0  =  g(Ab)  curve,  namely: 

A0  =  0,  Ab  =  +  Ab 

max 

A0  =  0,  Ab  =  -  Ab 

max 

Ab  =  0,  A0  =  +  A0 

max 

Ab  =  0,  A0  =  -  A0 

max 

A  first  approximation  to  the  A0  =  g(Ab)  curve  can  therefore  be 

obtained  by  simply  connecting  these  points  with  straight  lines. 

The  following  sketch  illustrates  this. 
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The   area    enclosed    by    this    first    approximation    is 

S    =   %[2(Ab         )A0   =    O][2b(A0  )Ab   =    o]    =   ^   S 

v      max'  v       max'  max 


which    leads    to   P    =   \  P 


max 


4d 


and  a  +(0)  =  -~~   [a  +(0)]  . 


(9) 


H 


o 


Equation  (9)  predicts  the  experimentally  observed  cross  section 

+ 
for  H   to  within  the  accuracy  of  the  experiment. 


61 


APPENDIX  A 


1.   Computer  Program  for  0"   (0)  and  CT   (9) 

H+  Ho 


This    computer    program   solves    the   following   set    of    equations 
for    classical    differential    scattering    in    the    center    of   mass 
coordinate   system. 


wher  e 


a     (o 

H 


db 


Sin9      d0 


Q    =   n     -    2bU 


1     _    V(r)        _    b_ 
•■,22 
%UVQ  r 


dr 


with 


r      =    largest  p 

o  r 


*.      *    i         V(r  )  b  _ 

ositive   root    of    1    -   — » — <-^    -   —^  -   0 

huv  r 

o 


11   =    reduced  mass    of    the    syst 


em 


V   =  velocity  of  incident  particle  in  the  laboratory 
system 
V(r)  =  scattering  potential 

In  this  case  the  scattering  potential  as  given  by  Byatt  L  3 J 

OiR         BR   VR 
for  Argon  converted  to  electron  volts  is  V(R)  =  -(Ae   +BeH  +C   ) 

K 

with  R  in  units  of  Bohr  radii  (0.5292A)  where 

A  =  411.4  &   =  -  1=676 

B  =-117=5  p  =  -  3.127 

C  =  195.9  y  =    -   9.625 
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An  outline  of  this  program  is  as  follows: 

1.  Read  in  the  desired  incident  energies,  mass  of  proton,  mass 
of  target  atom,  disassociat ion  energy,  and  potential  parameters. 

2.  Select  an  impact  parameter,  b.   From  an  initial  value,  b  is 


increased  in  units  of  .05  a   where  a   is  the  first  Bohr  radius. 

o         o 

3-    Solve    for    the   largest   positive   root   of 

.  2 


V(r  ) 
v  o' 


2        2 
/iV   /2    r 
o        o 


=  0  using  the  Newton -Raphs on  method. 


4-  Evaluate  the  integral 


2  2    J  2 

fLV  r  r 


using    the  Gauss -Mahler    quadrature   formula 


TT 

u  =  — 

r 
o 


(cos    ^)g(x1)  +  (cos   y^-)g(x2)  +  (cos   ^jOgC^) 


where   g(x . )    = 


1    - 


V(r    /x.) 
v    o      x' 

i  2 


x    o      i' 


11  3^  5^ 

with   x      =    cos   — -,    x      =    cos    — — ,    x      =    —    .       This   procedure  was    taken 

from    Bates,    Cook,    and  Smith   L  2  J . 


5.    From   0    =  tt    -    2bU,    evaluate    |— |   =    1 2b   +    2b  —  |.      —    i 


dU  i    dU 


obtained  from  the  analytic  expressions 

.  2 


U    =    f(r     ,b)    and    1    - 
v    o'     ' 


V(r     ) 
v    o' 

2  2 

11/2   V  r 

n  o  o 


db  '        db 


=   0    in    a    straight-forward 


dr 


manner    using    the   chain-rule   of   differentiating   —  =   —   +  r— —  -77—    . 

o 


i  db 


6.       Solve   for   a      (Q)    =    -^   |  —  | 
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Op. 

7.  Solve    for    the   orientation   probability,    P    =   — 0"       (0) 

tt£    E  H+ 

o 

8.  Solve    for    CT       (0)    =    P   Or       (9). 

H+  H+ 

2 


9°    Write   out    incident    energy,    b,    r     ,    9 ,    CT       (9),    P,    CT       (9) 

H+  HJ 


Computations  are  done  in  units  of  Bohr  radii,  but  CT   (9)  an 


d 


H2 


2 
CT       (9)    are    converted    to    cm      at    the    end   of    the   program.       The 

H2 
following    is    a    list    of    key    symbols    used    in    the    computer    program; 

DENER    =    Dfsassociatiar Energy ,     D,     in    units    of    eV. 
RHO(J)    =    Impact   parameter,    b,     in    units    of   Bohr    radii. 
EV(K)    =    Incident    ion    energy    in    units    of    eV . 

ET   =   %/l  EV(K>       definition 

H+ 

V(  R )  R 
V  =  — ' '  ^    convenient  form  of  the  potential  in 

this  program 
dV 

DV  ■  s 

DEG(J)  =  scattering  angle  in  degrees  converted  to 

laboratory  system 

SIGMA  =0"   (9) 
H+ 

cs  =  a     (0) 

H2 
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2 .   Calculation  of  Ion  Trajectory  and  Solid  Angle. 

As  discussed  in  Section  III  E,  the  scattered  ion  trajectories 
and  solid  angles  must  be  computed  numerically .   Since  much  of  the 
computer  program  that  accomplished  this  is  peculiar  to  the  Naval 
Postgraduate  School  IBM  36O-67  computer  installation,  only  an 
outline  of  the  program  is  presented  here. 

Recall  that  the  axial  component  of  the  magnetic  field  on  the 
magnetic  field  axis  was  measured  and  expressed  as  a  polynomial 


of  the  form 


13 

B(z,r=0)=IE   Cz 

z  \  n 


n-1 


Bush   L4J    shows    that    a    vector    potential    for    this    axially    symmetric, 
non-uniform   magnetic    field    can    be   written    using    this    expression 
for    B    ( z , r    =    0)    as; 


2n. 


UL 


Vr'2)    r    n=0 


2n+l      b       B    (z9r    =    0) 
r  7 v  ' 

(I)  ~ 


dz 


2n 


Ar(r,z)    =    0 


A    (r,z)    =    0 
z  v 

Using    this    vector    potential    in    the  Lorentz    force   equation 

F    =    qvxB=qv    x(V    x   A)  , 


it  follows  that 
,2 


F   =  m 
r 


F0  =  m 


F'  =  m 
z 


,(A4\ 


d  r 

dt2     ~\^> 

'dr 


A2 

d  2 


dt 


[d0  ^r  v" 

=  q 

dt      3r 

(1) 

2    i 

d  0 

r/dz/M   1  dr  a<r  vl 

j^2 

=  -  q 

(  dt  J\  dz    J  '    r    dz         dr 

dt     J 

\       /\  j     /                                    J 

=  q 

\    d0   ^' 
dt      $z  _ 

(3) 

(2) 


68 


d0 


^0 


Equation  (2)  can  be  reduced  to  — -  =  -  q  — —  .   Using  this, 

M        v  '  dt  ^  mr         " 

equations  (1)  and  (3)  can  be  combined  to  eliminate  the  time 
dependence  to 


dz2 


-  <§ >' 


A 


0 


dr 
dz 


dA 


A 


0 


0     dz 


••»•£-<■■ 


(4) 


This    is    a    second-order,    non-linear    differential    equation   whose 
solution    describes    the   trajectory    of    the   scattered    ion    in    the 
magnetic    field.       In    order    to    integrate    it    numerically,    it    is 
expressed   as    two    first-order,    simultaneous    differential    equations 
as    follows . 


Let         yn    =    r 


y^  =  X  = 


dr 
dz 


dy. 


dz 


dA 


0 


f       = 


0    dz 


*—    -   A    2 
2         A0 

q 


dA, 


f_   = 


l0  dr 


—Z   ~    Art2 
q  0 


Then    equation    (4)    may    be   equivalently   written    as 


dz 


=  y. 


(5) 


^2.     f      3 

dz    '         Iy2 


f  v      +  f  v      -    f 
2y2  ly2  2 


(6) 
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Equations  (5)  and  (6)  were  integrated  by  the  DHPGC  subroutine 
which  employs  a  Hamming's  modified  predictor -cor rector  method. 
The  solution  to  (5)  and  (6)  yields  the  trajectory  and  hence  Z 
and  0  .   Now  in  a  single  computer  run,  the  trajectory  equations 
were  solved  a  number  of  times  for  various  scattering  angles. 
Typically  this  would  be  for  0  from  35   to  55   in  1   increments. 
A  second  subroutine  was  then  called  on  to  generate  a  seventh 
order  least  squares  polynomial  fit  of  Q  =  f(Z  ) .   From  0   and 
the  detector  geometry,  AZ   was  computed.   The  solid  angle  could 
then  be  computed  by  the  expression 


°-  -®-in8i£§;i  K 


The  computer  output  then  tabulated  0,  Z  ,  AZ  ,  0  ,  and  d  CI    as 
well  as  a  plot  of  the  trajectory  (see  Fig°  12)  if  desired. 
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